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Abstract. This work has a purpose to collect selected facts about the com- 
pletely monotone (CM) functions that can be found in books and papers devoted 
to different areas of mathematics. We opted for lesser known ones, and for those 
which may help determining whether or not a given function is completely mono- 
tone. In particular, we emphasize the role of representation of a CM function as 
the Laplace transform of a measure, and we present and discuss a little known 
connection with log-convexity. Some of presented methods are illustrated by 
several examples involving Gamma and related functions. 
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1. Introduction 

A positive function defined on (0, +00) of the class C°°, such that the sequence of 
its derivatives alternates signs at every point, is called completely monotone (CM). 
A brief search in MathSciNet reveals total of 286 items that mention this class of 
functions in the title from 1932 till the end of the year 2011; 98 of them have been 
published since the beginning of 2006. 

This vintage topic was developed in 1920's/30's by S. Bernstein, F. Hausdorff and 
V. Widder, originally with relation to so called moment problem, cf. [3] [13j [TH [26j 
[27] . The much cited (but perhaps not that much read) Widder's book [28] contains 
a detailed account on properties of CM functions and their characterizations. The 
second volume of Feller's probability book [8 discusses CM functions through their 
relationship with infinitely divisible measures, which are fundamental in defining 
Levy processes. In past several decades, Levy processes have gained popularity in 
financial models, as well as in biology and physics; this is probably a reason for 
increased interest in CM functions, too. There are also other interesting topics 
in Probability and Statistics where CM functions play a role, see [TB] for one such 
topic. Aside from probability and measure theory, CM and related functions appear 
in the field of approximations of functions, as documented in the book [6] of 2007. 
Finally, they are naturally linked to various inequalities; several general inequalities 
for CM functions can be found in [17] , for a quite recent contribution in this area 
see [2]. 

This text has a purpose to collect well known facts about the CM functions, 
together with some less known ones, which may help determining whether or not 
a given function is completely monotone. In that sense, this work can be thought 
of as being an extension and supplement to another paper in the same spirit - |24] 
by Miller and Samko. In particular, we emphasize the role of representation of 
a CM function as a Laplace transform of a measure, and we present and discuss 
a little known (and even less being used) connection between CM function and 
log-convexity. Some of methods discussed in sections [2]-[5] are illustrated by several 
examples involving Gamma and related functions in Section [6] References and 
examples reflect author's preferences, and are by no means complete; the same can 
be said for the selection of topics that are discussed in this work. 
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2. Representations of completely monotone functions 

We start with a classical definition of CM functions, and we present two pos- 
sible representations in terms of integral transforms of measures and alternative 
representations for Steltjes transforms and CM probability densities. 

2.1 Definition. A function / defined on (0, +oo) is completely monotone if it 
has derivatives of all orders and 

(1) (-l) k f (k) (t) > 0, t G (0,+oo),fc = 0,1,2,... □ 

In particular, this implies that each CM on (0, +00) is positive, decreasing and 
convex, with concave first derivative. 

2.2 Limit properties. By (JXJ) , there exist limits of f( k '(x) as x — > for any 

k > 0; if those limits are finite, then / can be extended to [0, +00) and (Q} will also 
hold for x = (with strict inequality for all k). Limits at zero need not be finite, 
as in f(x) = -, for example. 

Clearly, lim a; _j. +00 f^(x) = for all k > 1. The limit of f(x) at +00 must be 
finite, and if it is non-zero, then it has to be positive (for example, f(x) = l+e~ x ). 

2.3 Lemma. The function f is CM if and only if [28] 

(2) /(*)=/ e-^d/x(t), 

•/[0,+oc) 

where fi(t) is a positive measure on Borel sets of [0,+oo) (that is, n(B) > for 
every Borel set B G R+ ) and the integral converges for < x < +00. 

In other words, completely monotone functions are real one-side Laplace trans- 
forms of a positive measure on [0,+oo). If the measure /j, has an atom at t = 0, 
then lim^-^+oo f(x) > 0. The measure [i is a probability measure if and only if 
lim :E ^o + /(0) = 1 (by monotone convergence theorem). 

The Lebesgue integral in can be expressed as a Lebesgue-Stieltjes integral 

(3) /(*)=/ e- xt dg(t), 

J[0,+oo) 

where g(t) = /i([0,£]) is the distribution function of /i, with <?(0_) = 0. For a 
positive measure the function g is non-decreasing, and by change of variables 
t = — log s we get 

2.4 Lemma. The function f is completely monotone on (0, +00) if and only if 

(4) f(x) = / s* dh(s), 

J[o,i] 

where h(s) = — g{— logs) is a non- decreasing function. 

If / is a CM which is the Laplace transform of a measure /i, as in ([5]), we 
write / = £( d[i) or f(x) = C{ dfi(t)). Similarly, the relation ([3]) between / and a 
distribution function g, can be denoted as / = L{ dg). If /1 has a density h with 
respect to Lebesgue measure, we write f(x) — C(h(t) dt) or only / = C{h). It 
follows from inversion formulae that each CM f determines one positive measure 
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fi via relation / = C{ d/i) and it is of interest in many applications to find that 
measure. 

2.5 Remark. Since measures are determined by their Laplace transforms, if 
f = £( d/x), then / is CM if and only if /i is a positive measure. If there exists a 
continuous density h of /i, then / is CM if and only if h(t) > for all t > O.D 

Let us now observe a subclass of CM functions which contains all functions / 
that can be represented as Stieltjes transform of some positive measure /z, that is, 

'[0,+Do) X + S 

It is easy to verify that each function of the form ([5]) with a positive measure /i 
is CM, hence / = C{v), where v is a positive measure. To find v, we start with 

— !— = { e -(*+s)« du, 

x + S J[0,+<x>) 

and, after a change of order of integration we arrive at the following result. 

2.6 Lemma. The Stieltjes transform of a positive measure fi as defined by |5P 
can be represented as a Laplace transform 



(5) /(*) 




J[0,+°o) \J[0,+oo) J 

That is, / = C{v), where the measure v is absolutely continuous with respect to 
Lebesgue measure, with a density £( d/i). 

Stieltjes transforms / have the property that —/is reciprocally convex (in ter- 
minology introduced in [21] , a function g[x) is reciprocally convex if it is defined 
for x > and concave there, whereas g(l/x) is convex). As proved in |21j . each 
reciprocally convex function generates an increasing sequence of quasi- arithmetic 
means, and hence CM functions that are also Stieltjes transforms are interesting as 
a tool for generating means. 

2.7 Completely monotone probability densities. Let / be a probability 
density with respect to Lebesgue measure on [0, +oo), that is, 

/ f(x) dx = 1 and f(x) > for all x > 0. 
Jo 

Then / is a CM function if and only if @ holds, which, after integration with 
respect to x £ (0, +oo) gives (via Fubini theorem for / > 0) 

r+oo i 

i-jf - fMt y 

Defining a new probability measure v by v{B) = f„ - dfj,(t), we have that 
(6) f(x) = f te- xt du(t) = [ te- xt dG(t), 

J[0,+oo) J[0,+oo) 

where G is the distribution function for v. The function x te~ xt is the density 
of exponential distribution Exp (t). Therefore, a density / of a probability measure 
on (0, +oo) is a CM function if and only if it is a mixture of exponential densities. 
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Note that (j6]) can be written as f(x) = E(Te xT ), where T is a random variable 
with distribution function G; by letting S = 1/T we find that 

(7) f(x)=E (L*) = ( -e^'dH(s), 

where H is the distribution function of S. The latter form is taken as a definition 
of what is meant by a CM density in [17l 18.B.5]; this is more natural than ([6]) 
because the mixing measure H is defined on values of expectations (s) of exponential 
distributions in the mixture, rather then on their reciprocal values as in ([5]). 



3. Further properties and connection with infinitely divisible 

measures 

Starting from the mentioned representations of CM functions, an interesting 
criterion for equality of two CM functions is derived in [7] : 

3.1 Lemma. If f and g are CM functions and if f(x n ) — g(x n ) for a positive 
sequence {x n } such that the series l/x n diverges, then f(x) = g(x) for all 
x>0. 

As a corollary to Lemma 13.11 we can see that if CM functions / and g agree 
in any subinterval of (0, +oo), then f(x) = g(x) for all x > 0. A converse result, 
which is also proved in [7] is more surprising: If / is CM and if the series J2n V x n 
converges, then there exists another CM function g ^ /, such that f(x n ) ^ g{x n ) 
for all n. 

3.2 Convolution and infinitely divisible measures. Given measures /i and 
v on [0, +oo) and their distribution functions g^ and g Vl we define the convolution 
/i * v as a measure with the distribution function defined by 

(8) 9n*v(t)= 9n(t - u) dg v (u) = / g u {t - v) dg^v) 

J[0,t] J[0,t] 

To show equality of integrals above, we use the formula for integration by parts in 
Lebesgue-Stieltjes integral (see [TB] or [3]) and note that the function u h-> g^t — u) 
is continuous from the left, while u i-> g v (u) is continuous from the right, hence the 
additional term due to discontinuities in the integration by parts formula equals 
zero, that is, 

/ 9^{t-u)dg„{u) = - \ g v {u) dg^t - u) 
J[o,t] J[o,t] 

and then we apply change of variables in the last integral, u = t — v. 

Repeated convolution is defined by induction, using associativity. In particular, 
the nth convolution power of a measure [A, denoted by fj, n * is defined by n — 1 
repeated convolutions [i * [i * ■ ■ ■ * [i. 

A measure fi is called infinitely divisible (ID) if for every natural number n there 
exists a measure [i n such that \x = /i™*. 

In the next two lemmas we collect some basic properties of CM functions. For a 
collection of other properties we refer to [21] . 

3.3 Lemma. If f and g are CM functions with f — C( dfi) and g = £( Av), then 
for a > 0, 

af = £(d(ari), f + g = C{d(p + v)), fg = C(d(jM*v)). 
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Therefore, if f,g are CM then af + bg (a, b > 0) and fg are also CM. 

Proof. First two properties follow from the definition of Laplace transform. 
The third property for arbitrary positive measures is proved in 8, p. 434]. 

3.4 Lemma, (i) If g' is CM, then the function x H> f(x) — e~ 9 ^ is CM. 
(ii) If log f is CM, then f is CM (the converse is not true). 
(Hi) If f is CM and g is a positive function with a CM derivative, then x > 
f{g{x)) is CM. 

Proof. To prove (i), let h(x) = er 9 ^ and note that h > and h! = -g'h < 0. 
Then by induction, using Leibniz chain rule, it follows that (— l) n h^ > 0. In 
particular, if log / is CM, then (—log/)' is also CM, and (ii) follows from (i) with 
g = — log/. The function x n- e~ x is a CM function but its logarithm is not the 
one, so the converse does not hold. For (iii), we note that / = £( d/u) for some 
positive measure /i, hence 
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(9) —f(g(x)) = -g'(x)J o e'^Hd^t) 

By part (i), the function x i-> e -9 ^* is CM for every t > 0, and so the function x i— > 

g' (x)e~ 9 ^ 1 is also CM as a product of two CM functions. Then from representation 

(J9j) it follows that the first derivative of —f(g(x)) is CM, which together with 

positivity of / yields the desired assertion. □ 

Note that if we can find measures v and v t in representations g'(x) = C{ Av) and 
e -g(x)t = dt , t ) ; then from jjgj we find that 



+ 00 Z' + OO 

'/ 
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(10) A.f(g( X ))=- I // ,-»•' (1(7/ *//,)(« ) (!//(/). 



It turns out that CM functions / of the form as in (i) of Lemma 13.41 are Laplace 
transforms of ID measures. If /(0) = 1, the associated measure is a probability 
measure, which is the case that is of interest in applications. Proofs of statements 
of the next lemma can be found in [8] . 

3.5 Lemma, (i) A function f is the Laplace transform of an ID probability 
measure if and only if 



(11) fix) 



where g is a positive function with a CM derivative and g(0) = 0. Equivalently, f 
is the Laplace transform of an id positive measure if and only if f(x) > for all 
x > 0, and the function x i— > — log/(ir) has a CM derivative. This measure is a 
probability measure if and only if /(0+) = 1. 

(ii) A function f is the the Laplace transform of an ID probability measure if 
and only if 

(12) -log/(x)= / d M (i), 

Jo 1 

where fi is a positive measure such that 



(13) / - d/i(t) < +oo. 
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3.6 Remarks. 1° If log / is CM, then — log / has a CM derivative and by Lemma 
I3.5( i). / = C( d/i), where /z is an ID positive measure. By CM property of log/, 
we have that log / = C(dv), where v is some other positive measure. Note that 
positivity of v implies that log/(0) > 0, that is, ^([0,+oo)) = /(0) > 1 and so, [i 
can not be a probability measure. 

2°. Non-negative functions with a CM first derivative have a special name - 
Bernstein functions; Lemmas 13.41 and 13.51 explain their role in probability theory; 
more about this class of functions can be found in |25j . 

4. M AJORIZATION, CONVEXITY AND LOGARITHMIC CONVEXITY 

A good source for studying all three topics that are very much interlaced, is the 
book [17] . In this short digest we include only necessary definitions and results that 
one can need for understanding a connection with CM functions. 

4.1 Majorization and Schur-convexity. For a vector x g R" define Xui to 
be the ith largest coordinate of x, so that 

X[l] > X[2] > ■ ■ • > X[n]- 

We say that x is majorized by y in notation x -< y if 

k k n n 

for fc= l,2,...,n- 1 and ^ = ^ 

i—1 i—1 i—1 i—1 

For example, (1,1,1) -< (2,1,0). Clearly, majorization is invariant to permuta- 
tions of coordinates of vectors. 

A function / which is defined on a symmetric set S C R" (S is symmetric 
if x G S implies that y € S where y is any vector obtained by permuting the 
coordinates of x) is called Schur-convex if for any x, y G S, 

(14) x -< y => /(x) < /(y). 

The following result, due to A. M. Fink [9, reveals an interesting relationship 
between concepts of Schur-convexity and complete monotonicity. 

4.2 Lemma. For a CM function f and a non-negative integer vector m = 
(mi, m2, • • • , rrid) of a dimension d > 1 . let 

u x (m) = (-l) mi /( mi >(x)(-l) m2 / (m2) '(as) ■ ■ ■ {-l) md f {md) {x). 

Then ^(m) is a Schur-convex function on m for every x > and d > 1. 

An important corollary of 14.21 is with d = 2, taking m = (1, 1) and n = (2, 0). 
Clearly, m -< n and from the above definition of Schur-convexity we get that 
u x (l, 1) < 1^(0,2), that is, (f'{x)) 2 < f(x)f"(x), which is, knowing that f(x) > 0, 
equivalent to (log f(x))" > 0. We formulate this result as a separate lemma. 

4.3 Lemma. Any CM function f is log-convex, i.e., the function \ogf(x) is 
convex. 

A converse does not hold, for example the Gamma function restricted to (0, +oo) 
is log-convex, but it is not CM . However, the fact that each CM function is also 
log-convex, helps us to search for possible candidates for complete monotonicity 
only among functions that are log-convex. In addition, there is a very rich theory 
that produces inequalities using convexity or Schur-convexity, and we can use it for 
CM functions. 
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Log-convexity of CM functions is equivalent to decreasing of the ratio f'(x)/f(x), 
and (arguing that /( 2fc ) and — j( 2fc + 1 ) are CM) this implies 

4.4 Corollary. If f is a CM function, then the ratio 

/(*+') (s) 
/(*)(*) 

is decreasing for every integers k,j. 

In the next lemma we give two consequences of convexity and log-convexity of 
CM functions. Similar inequalities for CM functions can be found in [16] , but with 
more involved proofs. 

4.5 Lemma. If f is completely monotone, then 



(15) f(x) + f(y)< f(x-e) + f(y + e) < /(0) + f(x + y), 

(16) f(x)f(y)< f(x-e)f(y + e) <f(0)f(x + y) 

where < e < x < y, assuming that /(0) is defined as f{x + ) (as in \2.'£\, finite or 
not). 

Proof. If f is a convex function, then the divided difference 

&<P,e{ X ) = 

is increasing with x, hence in the present setup, A.t^ £ (x) < Af^ e (y+e) and Ay i2 ,_ e (a;— 
e) < Af x ^ e (x + y), which proves (fT5|) . The same proof holds for (fTTj)) . but with 
log / in place of /. □ 

Let us note that under assumptions of Lemma 14.51 (x,y) -< (x — e,y + e) -< 
(0, x + y), and so we have just proved that the functions (x, y) i-> f(x) + f(y) and 
(x,y) i — y f(x)f(y) are Schur-convex on R + x R + . More generally, for any / being 
CM , the functions of n variables 

n n 

(17) and Y[f( Xi ) 

i=l i=l 

are Schur-convex on R™ . For a proof of this statement see [IT] . 

Finally, the fact that /' is concave (i.e, /'" < 0) is equivalent to each of three 
inequalities in the next lemma [TT)1 120) . 

4.6 Lemma. For a CM function f, it holds 

(18) < < / — — , for all x 7 y>0 



, 1fn f(y) ~ f(x) ^ f(y - g) - f(x + g) n ^ ^ 

(19) < , /or < a; < y and < e < ^^5— 

y — x y — x — 2e z 

5. Inversion formulae 

It is sometimes easier to find a measure /i that corresponds to function / via 
Laplace transform in ([3]) then to show that / is CM by verifying the definition; 
in view of applications, it is definitely useful and desirable to know the associated 
measure. In many cases we can use properties of Laplace transform and the tables 
that can be found in textbooks. In many applications the Laplace transform is 
not limited to real argument, and it is more common to define f(z) by ([3]), where 
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complex argument z belongs to some half space > a, for some positive a. We 
may use the power of complex Laplace transform calculus applied to real function 
of real argument, due to well known properties of regular functions. 

Due to similarity between Fourier transform, complex Laplace transform and 
real Laplace transform, we may use inversion formulae for all three mentioned 
classes, whenever it is appropriate. In probability theory, for a random variable 

the function x i— > E e txZ (which corresponds to Fourier transform, except the 
sign in the exponent) is called the characteristic function, whereas the real Laplace 
transform (mind the sign!) x i— > Ee lZ is called the moment generating function. 
There are several formulas that can be found in textbooks, but we will mention here 
only a not widely known inversion theorem that enables finding a finite measure fj, 
defined on Borcl sets of R, provided that we know its characteristic function 

r+oc 

(20) <p(x) = / e ltx dF(t), 



where Fit) = fi{(— oo, t]}. The following result (given here in a slightly generalized 
version) is due to Gil-Pelaez [12] . 

5.1 Lemma. For <p and F as in &20\) , with f(0) being finite, we have that, for 
all t e R, 

(2D m±iM = m _ i r R f^m\ dx . 



Note that the underlying measure here need not necessarily be restricted to the 
positive part of the real axis. As an example of how (|2"Tj) can be used to determine a 
measure \x such that / = C(fi), consider a simple case f{x) — e~ x , where we already 
know that the measure is Dirac at t = 1. Supposing that we wish to use (|2"Tj) to 
derive this, note that if / is the Laplace transform of /z, then its characteristic 
function is (f(x) = f(—ix) = e lx , and (|2"Tj) yields (assuming that t is a point of 
continuity of F) 

2 7r ,/n x 



Knowing that 



+ CXD 



sin(ax) ir 

da; = 77 s S n a i 



o 



x 2 

we find that F(t) = for t < 1 and F(t) = 1 for t > 1, hence (by right-continuity 
and non-decreasing of F), the corresponding measure is indeed a Dirac measure 
at t = 1. 

For other formulas and methods, including numerical evaluation of inverse, see 
[5]. In the next lemma we complement some examples from 24 by effectively 
finding the corresponding measure. 

5.2 Lemma. We have the following representations: 

(23) e- ax = £(dS a (t)), 

where 5 a is the probability measure with unit mass (Dirac measure) at a > ; 

(24) 1 c(e- bt ' a ^-\, a,b,c>0,a 2 + b 2 >0; 
(ax + b) c \ a c l (c) J 
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(25) log (a+ =£(d;u(i)) o>l,6>0, 

where the measure /i is determined by its distribution function 

f x 1 - e~ hs l a 

H([0,t]) =loga+ / ds; 

Jo s 

(26) ^i±^=£(^ (i )), 

x 

where (see [T| p. 56]) E\ is exponential integral 

/■ + °° dn 
Ji u 



(27) 



where I\ is a modified Bessel function as defined in pQ. 

Proof. The relation (|2"5|) is obvious, and (IM1) is a consequence of standard 
rules for (complex) Laplace transform: 

* C_1 \ _ 1 ruc-U, ,/ s 1 r(c) _ 1 



\ a c L (c) / a c l (c) 



o c r(c) (x-6/a) c (air + &) c 
To prove (|2"5")) . denote its left side by /, and observe that, by 

/'(*) = -1- - i = £ (e"^ - l) . 
ax + b x \ ) 

Now we use the rule 

C ( 9 M]=l C(g(t))[y]dy 



t J x 
to conclude that 

f(x) = log a - J f'(y) dy = £(log a dd ) - C ( J 



which yields (|25|) . To prove (|27jl . we note that 



a k „ ( a k t k 1 



= C 



which tells us that 



Now we observe that 



and (J27]) follows. 



klx k \M(k- 1)1 J' 



1+£ £ 



' (fc - l)!jfe! 



^ a k t k - 1 _ ahi^M) 
^ (k - l)\k\ ~ Vol ' 
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The simplest way to prove (|26j) would be to perform an integration on the right 
hand side and show that it yields the left side. However, in order to show the 
derivation, we start with the observation that 

f( X ) :=^±^=F(l,l,2;-,), 

x 

where F(a, b, c; •) = 2^1(0, b, c; •) is Gauss' hypergeometric function; hence there is 
the following integral representation pQ : 

r 1 ds 



m = 

Now we use (f2"5|) to find that 

1 1 



1 + sx 



+00 

e -xt e -t/s dt 



1 + SX s 

and, exchanging the order of integration, we find that 



+00 / ,1 d 



/(»)= e-* e-*/'^ dt. 



Finally, a change of variables 1/s = it in the inner integral shows that it is equal to 
Ei(t), and the formula is proved. 



6. Some examples related to the Gamma function 

Functions related to the Gamma function are good candidates to be CM , and 
there is a plenty of such results in literature. The function g(x) — logr(a;) is a 
unique convex solution of the Krull's functional equation 

(28) g(x + 1) - g(x) = f(x), x > 0, 

with f(x) = log a; and with g(l) — 0. The same equation, but with f(x) — 
(log£)(™ +1 ), n = 0, 1,2, . . . has for its solutions functions $W(x) = (logr(a;))(™+ 1 ). 
Although logx is not CM, all its derivatives are monotone functions, which auto- 
matically implies the same property for ty( n '(x),n > 2 and alike functions via the 
following result (see [22]). 

6.1 Lemma. Suppose that x 1— > f(x) is a function of the class C°°(0, +00) with 
all derivatives being monotone functions, with f'(x) —¥ as x — > +00. Then there 
is a unique (up to an additive constant) solution g of H28\) in the class C°° , with 



(29) g'(x)= lim [f(x + n) - V f'(x + k) 

and 



n— ^+00 



-\- 00 



(30) gW(x)=-J2f U) (x + k) (j>2). 



k=0 



From ((29j) and ([30]) it follows that, if ±f is CM (or if only ±f" is such), then 
=F</' is a CM, while ±g and ±g' need not be CM. Our first example is formulated 
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in the form of a lemma, and its proof provides a pattern that can be used in many 
similar cases. 

6.2 Lemma. The function 

W(x) = ~{\ogT{x) -(x- l)loga;)" = - + — - 



x x 2 



has the following integral representation 



(31) W(x) = ^ [1 + t- Y— t ) e- xt dt 

and it is a CM function . 

Proof. The integral representation follows from 

(32) c({-l) n + 1 T ^] =*W(aO, 71=1,2,..., 
and 

(33) £( t a ) = r^l) i fl> _ L 

The CM property follows from positivity of the function under integral sign, which 
is equivalent to the inequality e* > 1 + 1 for t > 0. 

Remark. The function g(a;) = logr(x) — (x — l)loga; satisfies the functional 
equation (|2"8)l with f(x) — log f jqpr) I ^ can be easily checked that /" is CM, 
hence from Lemma |6. II we can conclude without any additional work that — W as 
defined above is CM . 

6.3 Example. For a > and x > 0, fei 



( X ~ 5) l0g ' T ~ 12*'^ + a) + x ~ \ log ( 2?r )- 



G„(a:)= logics)- he 
T/ie following representation holds: 

(34) GaM ^- t -^(^,^-" e -., d , 

The function x G a (x) is CM if and only if a > 1/2 and the function x 1— » — G a {x) 
is CM if and only if a = . 
Proof. Starting with 

G2(z) = *'(*) - + a) - I + i^/?, 

it easy to show (in a similar way as in Lemma l6.2p that 

(35) . r t - 2+(2 :r'- (ta/6)e '" e- * 



Further, we have that 



2(1 -e-*) 



lim G (x) = lim G„(x) = 0, 

X— >-\-OQ X— + 



and 

G a (x) = / / G«(«)d«d«, 
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hence ([34[) holds. The complete monotonicity is related to the sign of the function 

t 3 

(36) h a (t) = t- 2+ (2 + t)e -i - — e~ at . 

The function G a is CM if and only if h a (t) > for all t > 0. From ((551) we see that 
this is equivalent to 

(37) a > los6 + log((2 + t)e-* +f -2)-31og f ;= ^ 

Using standard methods, we can find that u is a decreasing function, hence 

u(t) < lim u(t) = -, 
t— rj+ ^ 

and so, (l37l) holds if and only if a > 1/2. 

Further, — G a is CM if and only if h a (t) < for all t > 0, which is equivalent to 

(38) a < u(t), 

where w(f) is defined in ((37]). Since w is decreasing, we have that 

u(t) > lim u(t) = 0, 

t— f+oo 

and so, ([55)1 holds if and only if a < 0, that is, a = 0. 
Remark. Let 

F a (x) = logr(x) - he - logs- ^*'(x + a), a > 0,x > 

This function is studied in [181 Theorem 1], where it is shown that x i— >• Fq(x) 
is concave on x > and that x H> F a (a;) is convex on x > for a > i. Since 
F a {x)" = G"{x) where G a is defined as above, this example gives much stronger 
statement. 

6.4 Example. For b >0 and c > 0, let 

(39) / 6 , c (x) = eXT ^ b \ x > 0. 
The function 

(40) ^b,c(a:) = log/6, c (x) = x + logT(x + b)- (x + c) logx 

is CM if and only if b > | + — and c = — \ and then it has the representation 

(«) + * >a 

Proof. By expanding logr(a; + b) in (14T)[) by means of Stirling's formula [TJ 
p. 258], it follows that, for 6 = b — c ^ |, 

lim <£ b , c (x) = ( 5- \ ] • (+oo), 

so ¥>b, c is not a CM function (see I2.2[) . Let 5 = 1/2 and let 

i — 

2 



G 6 (a;) := V bb _i(x) =x + logT(.x + 6) - \ x + b- - )\ogx. 



Further, we find without difficulties that 

(42) lim G b (x) = lim G' b (x) = 
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and that 



G'l{x) = 9\x + b) - 




In the same way as shown in Lemma \Q.2l we find that G'l{x) = C(hb(t) dt), where 

(43) hb ® = T^rt +t ( b -l)- 1 - 

By standard methods we find that 

(44) h(t)= (?-- b - + ±^t>+o{t>) (t -►<)), 

so the Laplace transform Gj,(x) of the function t h-> g(t)/t 2 exists for all x > and 
applying Fubini theorem as in Example 16.31 and using (|42[) we find that 

G b (x) = [ °° / +0 ° G» d U dv= f ^e- tx dt, 

Jx Jv J[0. + oa) * 

which is the representation (|3T]) . Then Gf, will be CM if and only if if hit) > for 
each i > (see Remark |23|). By (j44|) we have that /i(0) < for b € 62), where 
&1.2 = I ± -7=; further, c = b — 1/2 > gives 6 > 1/2, so only 6 > 62 remains as a 

possibility. It is straightforward to check that 8 ^ > for all i > 0, so it suffices 
to show that h^it) > for t > 0, which can be done along the lines of |10) . 

Remark. Complete monotonicity of functions f b c and c for various values 
of parameters was discussed in [TD] and [TT]. Let us remark that, by Lemma \3. 41 
the function ft, >c is CM whenever (/Jfc.c is the one. □ 

Let us mention that the Barnes function G(a;) satisfies the relation 

log<3(af + l) -logG(x) =logr(^), a;>0, 

which is (J28]) with g = logT. Here also the function x ^ (logG(a;))" = 2^f'{x)+(x- 
l)*S?"(x) is CM. More details about the properties of the G-function as a solution 
of Krull's equation can be found in [23] . 
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